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We present a theoretical analysis of a microwave spectroscopy experiment on a charge qubit 
defined by a donor pair in silicon, for which we calculate Hamiltonian parameters using the 
effective-mass theory of shallow donors. We solve the master equation of the driven system in a 
dissipative environment to predict experimental outcomes. We describe how to calculate physical 
parameters of the system from such experimental results, including the dephasing time, T2, and 
the ratio of the resonant Rabi frequency to the relaxation rate. Finally we calculate probability 
distributions for experimentally relevant system parameters for a particular fabrication regime. 

I. INTRODUCTION 

The strong motivation supplied by the possibility of quantum information processing has recently fuelled rapid 
progress in the experimental control of mesoscopic quantum systems. Of particular interest in the solid-state are 
superconducting devices, and coupled quantum dots, which show promise as candidates for qubits in a scalable 
quantum device. In order to realise this potential it is necessary that these systems be coherently controlled with 
extraordinary precision, and much progress has been made toward this goal in superconducting systemsi^, as well 
as coupled quantum dots in gallium arsenide 4 , and phosphorus doped silicon^. 

In terms of the experimental demonstrations of coherent control one can take two approaches, the observation of 
Ramsey fringes in a pulsed experimenliiiiMi^, where the coherent oscillations are between non-diagonal eigenstates of 
the system Hamiltonian. The system must be manipulated on a timescale that is fast compared to the period of the 
oscillations, which is inversely proportional energy difference of the eigenstates. Alternatively one could observe the 
Rabi oscillations in a driven system^, where the frequency of the oscillations is proportional to the intensity of the 
driving field. 

In many cases a precursor to these achievements is the somewhat simpler procedure of performing a spectroscopic 
measurement of the devic o 6 i 7 i 8 . This is achieved by driving the system of interest at particular frequency, and varying 
parameters to bring the system into resonance. When in resonance, the system will undergo Rabi oscillations, however 
a spectroscopic experiment of this nature does not resolve these oscillations, instead measurements are made on a 
timescale that is greater than the Rabi frequency, giving a time averaged result. Useful results can be obtained even 
in the regime in which the frequency of the oscillation is slow compared with relevant decoherence times. Clearly 
a measurement of this type places far fewer demands on the system than if the fringes are resolved, however much 
useful information can be gained about the system, such as the dephasing time T2, some of the system energy scales, 
and the ratio of the resonant Rabi frequency to the relaxation rate. 

One of the most promising quantum computer architectures, particularly from the all- important scalability perspec- 
tive, is the phosphorus-doped silicon (Si:P) system where the shallow P-bound electronic donor states are electrically 
manipulated by external gates for the quantum control of the qubits. In spite of impressive recent growth, fabrica- 
tion, and theoretical developments in the Si:P quantum computer architecture, the key ingredient of an experimental 
demonstration of coherent quantum dynamics at the single qubit level is still lacking. In this paper we consider a 
specific experimental set-up, namely the externally driven microwave spectroscopy on a P^" charge qubitS. (i.e. the 
singly ionised coupled P-P effective 'molecular' system in Si with two nearby ion implanted P atoms sharing one elec- 
tron), giving a detailed quantitative theoretical analysis to establish the feasibility of this coherent single qubit control 
experiment in currently available Si:P devices. We believe that the qubit coherent control experiment proposed and 
analysed in this work may very well be the simplest single qubit measurement one could envision in the Si:P quantum 
computer architecture, and as such, our detailed quantitative analysis should motivate serious experimental efforts 
using microwave spectroscopy. The experiment being analysed in this work is a necessary prerequisite for further 
experimental progress in coherent qubit control and manipulation in the Si:P quantum computer architecture. Our 
detailed theory provides the quantitative constraints on the experimental system parameters required for developing 
the coherent control of the Si:P system, and in addition, establishes how one can extract important decoherence 
parameters (e.g. T2) from the experimental data. 

The energetics of this system, including its interaction with phonon modes in the lattice^, and with a DC electric 
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FIG. 1: (Colour online) On the left is a schematic of the experimental setup showing the location of the bias gates relative 
to the phosphorus donors. On the right is a diagram showing orientation of the two donors, and the direction of the applied 
fields, relative to the silicon crystal lattice. 

fieldii, have been studied by previous authors, while a theoretical investigation of the microwave spectroscopy of a 
coupled quantum dot in a GaAs system has been carried out by Barrett et afi£, and much work has been done on the 
theory of the dynamics of driven systems coupled to dissipative environments has been well studied in the context of 
flux qubitSfi&il 

In this article we begin by reviewing the in silicon system in the presence of an electric field, and outline the 
calculation of the Hamiltonian matrix elements for use in the calculation of the system dynamics. In sec lllll we derive 
a master equation for the driven system, in the presence of three decoherence channels, following the the approach 
of Barrett. This master equation is analytically solved for the steady-state solutions, valid on a time-scale that is 
long compared to the relaxation time. To probe the dynamics on shorter time-scales requires, in general, numerical 
integration of the master equation. In sec II VI we show how it is possible to use the output of the measurement device 
to obtain the energy scales of the system, and to calculate the dephasing time T 2 , while in sec. we briefly discuss 
possibilities for the direct observation of coherent oscillations in this system. Finally, in sec lVll we discuss the statistical 
nature of the device fabrication, we calculate probability distributions for experimentally relevant parameters, and 
discuss likely experimental outcomes. 
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II. A TWO-DONOR CHARGE QUBIT IN SILICON 

We consider a system defined by the electronic state of a pair of phosphorus donors in silicon, which has been singly 
ionised such that the system has a single net positive charge^. The donor nuclei are separated by a vector R, and the 
in the absence of externally applied fields, the system Hamiltonian is 

H Q = H Si + V d (r - R/2) + V d (r + R/2). (1) 

Here Hsi is the Hamiltonian of an electron in the pure silicon lattice, which includes both a kinetic term and the 
effective potential due to the silicon lattice. Solutions of this silicon Hamiltonian are the Bloch functions of the pure 
silicon crystal. The Vd(v) terms give the impurity potential of the ionic donor cores, which we will treat as Coulombic, 
Vd(r) = —2/nr, with n = 11.9, the dielectric constant of silicon. 

The electronic ground state for a single donor, centred at a position R/2 is, in the effective-mass approximation, 
given by the Kohn-Luttinger wavefunctioni^ 

^(r) = ^F /1 (r-R/2K (i (r)e*-( r - R / 2 ). (2) 

In this equation the sum is over the 6 degenerate minima of the silicon conduction band, denoting the reciprocal 
lattice vector at the minima fi. The function (r) is the periodic part of the conduction band Bloch function, and 
-F)i(r) is a non-isotropic, hydrogen-like envelope function, which for the ±z minima takes the form 

-Jz*/al+(x*+y*)/al 

F± z = == , (3) 

J 67ra^a|| 

where an, a± are the non-isotropic effective masses associated with the conduction band minima. 

If R is sufficiently large, the low-lying energy states of the two-centre system, are approximated by even and 
odd superpositions of the single donor ground-states, v I'±(r) = 1/^/2(1 + S 2 )(-ip(r — R/2) ± ip(r + R/2)), where 
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FIG. 2: (Colour online) A plot of the symmetric anti-symmetric energy gap A s _ as as a function of donor separation along three 
high symmetry crystallographic axes. The points denote fee lattice sites, the lines are a guide to the eye. 
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S = J t/)*(r — R/2)^(r + R/2)dr. Interference between Bloch functions at the different conduction-band minima leads 
to a non-monotonic dependence of the symmetric-antisymmetric energy gap, A s _ as , on the magnitude of the donor 
separation, as well as on their orientation with respect to the silicon substrateafi. This is illustrated in Fig. [21 where 
we plot A s _ as as a function of donor separation, for donors separated along three high symmetry crystallographic 
axes. Note that for certain separations, the ground state is actually the antisymmetric state. 
In the basis of the single donor ground states , we can write 

H = h x a x + h z a z , (4) 

where the coefficients H a are determined by numerical integration for various vectors R. Note we have ignored the 
constant term, and defined our axes such that the a y component is zero. In practice we expect that the two donors 
are identical, in which case H z = 0, however we retain this term in what follows for generality. 

The addition of a uniform electric field in a direction along a [100] axis, gives rise to a potential that can be written 

as 

V = Efoo* + (5) 

where E is the strength of the applied field, and the coefficients £ x ,xi z are determined by numerical integration. In 
the case of radially symmetric single donor wavefunctions, the £ x term is identically zero, however this symmetry is 
broken by the silicon lattice. In practice, we find, £ z >> £ x . The magnitude this difference is dependent on R, but 
for R w 20nm, it is typically a factor of between 100-1000. 
To begin with, we consider the effect of a DC field 

H = H + V dc 

= (K + SdcCxK + (h z + E dc £ z )a z 

= ecxV 0CT \ (6) 

Here = \J (h x + E dc ^ x ) 2 + (h z + E^z) 2 , and 4> = arctan [(h x + Edc£x)/(h z + Ed c £,z)]- Thus, the above Hamilto- 
nian can be diagonalised via the transformation: 

H' = UiHU\ = 6cr z , (7) 

with U\ = exp(i</>/2er y ). The eigenvalues are given by A± = ±0, and the eigenstates, in the original charge localised 
basis, by: 

|+) = cos[cA/2]|L)- S m[cA/2]|ii}, 

|-) = sin[c4/2]|L) + cos[c6/2]|i?). (8) 

Thus by varying the strength of the field we can alter the energy splitting A = 20, between the two lowest energy 
states, as well as the degree of localisation of the eigenstates \± = I < L\± > | 2 = 1/2(1 ± cos(cA)). 
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We now consider the additional of an AC field of angular frequency to applied to the uniform DC field. The 
Hamiltonian of this driving field is 

14c = E ac cos(ujt)^ x a x + i z a z ), (9) 
which transforms to the energy eigenbasis of the time independent system as 

= E ac cos(ut){a x lx +a z lz ), (10) 

with -y x = £ x cos((f>) — £ 2 sin (</>), j z = £ z cos(0) + £ x sin(0). We now transform into an interaction picture, rotating 
at the same frequency as the AC field, the transformation is defined by Uiit) = Q lu / 2t<y . Evolution in this frame is 
generated by the Hamiltonian 

H' = (6 - uj/2)a z + E ac cos(cjt)(a x lx e-^ + a z lz ). (11) 

It is, at this point that it becomes convenient to make a rotating wave approximation (RWA), in which we ignore 
terms that are rotating with frequency u. This approximation is valid if u) » E ac ^f x , i?ac7z> an d reduces the above 
Hamiltonian to 

& = (e-Lu/2)a z + E ac ^a x . (12) 



III. CONSTRUCTION OF THE MASTER EQUATION 

To include the effects of a dissipative environment on the dynamics of the system, it is necessary to construct 
a master equation. We consider the effects of three channels of decoherence, namely; dephasing, relaxation and 
excitation, with different rates r z ,T_,r + respectively. We expect these channels to be sufficient to describe most 
physical noise sourcesA&ii including interaction with the measurement devic o 12 i 18 . Generally we would expect, T z to 
be the dominant term, with the relaxation term T_ < T z . The excitation term is included for generality, and allows 
the description of the effects of classical, or high temperature noise sources, we expect T + < T_. 

These decoherence channels are included in the master equation, in the energy eigenbasis, as Linbladian terms 
£{Li, p} = Liph\ — l/2(L\Lip + pL[Li). thus, in the rotating frame, the master equation is given by 

P'{t) = -i [(6 - w/2)a* + £^y 

+ T z C{a z > p'(t)}+T_£{a-,p J (t)} + T + £{a+,p J (t)}, (13) 

The steady state solution of this master equation can be found, which gives the asymptotic form for the components 
of the polarisation vector, defined by p = 1/2(1 + Xa x + Ya y + Za z ), in the rotating frame: 

, = -877 7 (r + -r_) 

1 ' i6(r + + r_)7 ? 2 + (2 ( ^ ac7:c ) 2 + (r + + r_)(r + + r_ + 4r z ))(r + + r_+4r 2 ) 

2r ? 7(r + -r_)(r+ + r_ + 4r 2 ) 



y'(oo) = 



i6(r+ + r_)r? 2 + (2{E aclx y + (r+ + r_)(r+ + r_ + 4r,))(r+ + r_ + 4r, 



y l{ , (r _ r , ler? 2 + (r+ + r_ + 4r z ) 2 

[ ! [ + " j i6r ? 2 (r + + r_) + (2(i? ac72; ) 2 + (r + + r_)(r + + r_ + 4r z ))(r + + r_ + 4r z )- 1 ' 

where r; = — uj/2 is the detuning of the microwave field. 

We wish to calculate the expected output from a detector which measures charge localisation, a z-measurement in 
the original basis, and therefore are interested in calculating Z-componcnt of the polarisation vector in the charge 
localised basis and the laboratory frame, the asymptotic form of which is given by 

Z(t) = (Z' cos(0) - (cos(c^)A' + sm(ujt)Y') sin(0)). (15) 

In the case that the measurement time of the detector is long compared to both the relaxation timescale, idct > 1 /r_ , 
and the timescale of these oscillations ^et > l/fi 5 then the system reaches the steady state equilibrium over the time 
of the measurement, and the the rapid oscillations of the rotating frame will average out. In this case the probability 
of measuring the electron to be localised at the position of the left donor, is given by Pl = (1 + cos(</>)Z'(oo))/2. 
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IV. APPLICATION 



We now turn our attention to analysis of simulated experimental results. In practice the P-P + charge qubit is, in the 
first instance, fabricated via ion implantation^, whereby 14keV phosphorus ions are implanted into a silicon crystal. 
Such a process is imprecise and there will be significant uncertainty in the position of the two phosphorus donors. As 
we have previously discussed, this will be manifest in an uncertainty in both the bare Hamiltonian of the two donor 
system, as well as in the potential arising from the applied electric fields. To illustrate how this uncertainty will affect 
experimental results, we have calculated simulated experimental traces, Pl as a function of applied DC field -Edc, for 
donor pairs with with slightly varying separations, using a v — 40GHz driving field of amplitude E ac = 10~ 2 MVm . 
The results are shown in Fig. |3J and show that although a slight variation of the donor separation has little effect on 
the position of the resonance peak, both the width and height of the peak are strongly affected. 

In addition to the system Hamiltonian parameters, an experimental trace will be effected by the strength of the 
decoherence channels, as shown in Fig.^J Here we see that both the height, and the width of the peaks are dependant 
of the rates of decoherence, whilst the position of the peak, is determined by the eigen energies of the system. It is, 
therefore, reasonable to expect that useful information about the system can be obtained from the analysis of this 
spectroscopic data. 

Given an experimental trace, we would like to determine useful spectroscopic information about the qubit, in 
particular we would like to determine decoherence times, as well as the energy scales of the system. Observation of 
the height of the baseline, taken at a point where this baseline is flat, allows the evaluation of the ratio r+/T_. In 
the absence of microwaves, the signal will be Pl = 1/2(1 + (T+ — r_) cos(</>) / (T+ +T-)), which allows the evaluation 
of cos(4>) over the range of E& c . 

In general, the strength of the applied electrostatic fields at the qubit, both AC and DC, may not be known for 
a given control gate bias. In what follows we will assume that both these quantities are linear functions of their 
respective voltage biases, E^ c — aG<± cl E ac = /3G ac where the Gdc; G ac are known experimental control parameters. 

By obtaining several traces taken at different driving frequencies, and observing the positions of the resonance peaks 
Ej c , it is possible to obtain values for as a function of Gd c - Combined with the knowledge of cos</>, this allows us 
to obtain values for a£ z and h x . 

Additional information can be gained from measuring the width of the resonance peak, the half-width at half 
maximum, measured from the flat baseline value, is given by the expression 



AG dc = — / ( 2 (^c7*) 2 + (r+ + r_)(r + + r_ + 4r z ))(r + + r_ + 4r z ) 
2a£z y r + + r_ 

Here we have assumed that cos(</>) = 1 at the position of the resonance, which is a good approximation in all the 
traces shown in this article. If this is not the case, a slightly more complicated expression can be derived for Eq 1161 
If traces can be evaluated for different values of the microwave power, this expression gives information about the 
decoherence rates, in particular extrapolation to zero microwave power yields AGd c = l/(2a£ z )(r + +r_+4r z ). Given 
that a£ z has already been determined, this allows direct evaluation of the dephasing time T2 = 4/(T + + T_ + 4T Z ). 
An expression for the height of the resonance peak, 

1 ( (r + -r_)(r + + r_ + 4r z ) 

~2^ 2(£ ac7;c ) 2 + (r + + r_)(r + + r_+4r z ) 

where we have again assumed cos(</>) = 1 at resonance, allows the evaluation of the ratio r + /T_. 




V. MEASUREMENT OF COHERENT OSCILLATIONS 



A somewhat more demanding experiment would be to attempt to resolve coherent oscillations of the system. 
This can be done with either Rabi oscillations 2 , or Ramsey oscillations^^. Rabi oscillations are resonant driven 
oscillations of the type described in this paper, where the Rabi frequency is given by E ac j x . In order to be able 
to observe Rabi oscillations, it is necessary that the oscillation period be long compared to the time resolution of 
the detector tjet < ^/(E ac Jx), or alternatively, that the state can be shelved for measurement. In either case, Rabi 
frequency must be fast compared to the dephasing time E ac j x > T z + T + + T-. We have calculated a distribution 
of Rabi frequencies, in units of the applied AC field strength, for a certain fabrication strategy that is described 
in the following section, and plotted the results in Fig HJl. From this we can see, that for an AC field strength 
of E ac = O.OlMVm -1 , the value used for all calculations in this paper, and one which is possibly at the high end 
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FIG. 3: (Colour online) Donor spectroscopy for several donor separations with a small deviation from the ideal separation of 
R = [40a, 0, 0], where a = 0.543nm is the lattice constant of silicon. The data shows that the peak heights and positions are not 
strongly dependant on these positional variations, however the width of the peak changes significantly. All data were calculated 
with a 40GHz driving field of amplitude E ac = 10~ 2 MV/m, and decoherence rates F z = 3GHz,r_ = 1GHz and F+ = 100MHz. 




FIG. 4: (Colour online) An illustration of the effects of various decoherence channels on the spectroscopic signal height. These 
channels affect both the peak height and width as described in the text. The lengths are in units of a = 0.543nm, the lattice 
constant of silicon, all data in this plot was taken for a donor separation R = [40a, 0,0], where a = 0.543nm is the lattice 
constant of silicon. The frequency of the driving field was 40GHz and its amplitude E ac = O.OlMVm" 1 . The decoherence rates 
are a) F z = 3GHz, T_ = 1GHz, T+ = 100MHz, b)r z = 3GHz, T_ = 10MHz, T+ = 1MHz, c)r z = 300MHz, T_ = 10MHz, 
r+ = 1MHz, c)r z = 300MHz, T_ = 100MHz, F+ = 1MHz. 



of what could be reasonable expected to be experimentally generated, the Rabi frequency is almost certain to be 
significantly less than 500MHz. This is likely to be slow compared to realistic dephasing rates, making resolution of 
Rabi oscillations unlikely. 

For most quantum systems, and this case is no different, it may be more practical to attempt to measure the 
Ramsey oscillations. These are the coherent oscillation between basis states that are not eigenstates of the system 
Hamiltonian. The advantage observing Ramsey oscillation over Rabi oscillation is that the frequency of the Ramsey 
oscillations is significantly higher than the typical Rabi frequency. The frequency of the Ramsey fringes is given by 
the energy difference between eigenstates, 20, and can be increased with the application of a DC field E^c, to at 
least 40GHz. Apart from the obvious condition 20 > T z + T + + T_, and assuming that the system can be shelved 
for measurement, the experimental requirements for resolution of these oscillations are that the measurement time is 
fast compared to the relaxation time t^ ct < l/r_, and that the system can be manipulated on a timescale that is fast 
compared to the period of oscillations. 
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In General, Ramsey oscillations are significantly easier to observe than the Rabi oscillation, and have been observed 
in both superconducting systems 1 ^, as well as coupled quantum dotsi, whilst Rabi oscillations have only been resolved 
in superconducting systems^. 

VI. UNCERTAINTIES DUE TO IMPRECISE FABRICATION 

We now turn our attention to uncertainties that may arise due to the statistical nature of the fabrication procedure 
for such a device. Atomically precise placement of phosphorus donors in beyond current technology, even the most 
reliable fabrication procedures lead to an uncertainty of at least several lattice sites in the final positions of the 
donors^. As always, there is a tradeoff here between precision and speed, and the most practical way to create such 
a device is via the implantation of a pair of single phosphorus ionsi^, which leads to significantly higher uncertainty 
in the donor positions. 

Interactions of the implanted ion with the host silicon atoms gives rise to a phenomenon known as ion straggle, 
which results in an uncertainty in the final position of the ion. Implantation of 14keV phosphorus ions in silicon 
results in a roughly Gaussian donor distribution, centred approximately 20nm below the silicon surface, and with a 
deviation a w Ylnm. In addition, the relatively large spot size of the ion beam means that the ion is equally likely to 
enter the silicon substrate at any point in the mask aperture, which is approximately 30nm in diameter. The result 
is a complicated distribution of donor pairs, with a uniform distribution angles between the donor separation and the 
applied fields, decreasing the effective field strength. An alternative method would be to implant a single ion into two 
separate apertures aligned with the field. In Fig [5] we show that the single aperture strategy is the best approach. 
The Rabi frequency decreases as ^ x cos(9) exp(— i?/a), where 9 is the angle between the donor separation and the 
applied electric field, and so it is far more important to keep the donor separation to a minimum. 

Using a typical distribution calculated for a 14keV ion implant in a single 30nm aperture, we have calculated the 
distribution of certain parameters of the system which affect the spectroscopic trace, using an ensemble of 16000 
points. These results are plotted in Fig The distribution of the magnitude of the donor separation, is shown 
in Fig EJa), it is peaked around R = 30nm, and has a variance of the order of lOnm. The distribution of the 
unperturbed energy gap A is plotted in Fig E[b) , although it is strongly peaked at low frequency we see that there is 
a finite probability of finding A > 40GHz. In this case there is no chance of observing a single photon resonance peak 
using 40GHz microwave radiation, and although two-photon resonances may be observed, they will be of diminished 
magnitude. We also note here that we have plotted the magnitude of the quantity A s _ as plotted in refi^, and so do 
not observe the negative values seen by these authors. In Fig EJc) shows the distribution of E* lc , the DC field required 
to bring the system into resonance, which exhibits quite a large spread. This is partly to do with the distribution 
of 9, the angle between the donor separation and the applied electric field, clearly E^ c oc 1/ cos(0), and so for large 
angles may be experimentally unachievable. Fig IHId) shows the distribution of the quantity cos(0) at resonance. 
This distribution is strongly peaked around 1, indicating that at resonance the eigenstates are likely to be the charge 
localised states, and the approximation made in deriving expression Eq 1161 is likely to be valid. Finally, Fig [ffi el 
shows the distribution of Rabi frequencies, in units of the applied AC field strength, as discussed earlier, for achievable 
field strengths this is likely to be much lower than the frequency of Ramsey oscillations. 

VII. CONSLUSION 

We have developed in this paper a detailed quantitative theory for the electric field tuned microwave spectroscopic 
coherent response of the electronic P2 charge qubits in silicon. Our realistic calculation, taking into account fabrication 
and ion implantation aspects of the sample preparation, shows convincingly that such an electric field tuned coherent 
microwave response measurement to be feasible in the Si:P2 quantum computer architecture. In addition, we establish 
that such a measurement would directly provide information on the characteristic decoherence and visibility times for 
Si:P charge qubits. We find that the "exchange-oscillation" type quantum interference problem^, arising from the 
six- valley degeneracy in the Si conduction band, does not lead to any specific difficulties in the interpretation of the 
microwave spectroscopy proposed herein. Our proposed microwave experiment along with the recent dc experiment 
proposed by Koiller et ali^, which complements our proposal, are essential pre-requisites for the eventual observation 
of Rabi and Ramsey coherent oscillations in the Si quantum computer architecture. 
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FIG. 5: (Colour online) Donor spectroscopy for several donor separations with significantly different separations. This data 
plot is intended to illustrate that increasing the angle between the donor separation and the direction of the applied fields, 
has less effect on the spectroscopic system than increasing the donor separation. All data were calculated with a 40GHz 
driving field of amplitude E ac = 10 _2 MV/m. The decoherence rates for the separations R — [40a, 0, 0] and R = [40a, 40a, 0] 
are F z = 3GHz,r_ = 1GHz and T+ = 100MHz, whilst for R = [110a, 0,0] we have used F z = 300MHz,r_ = 100MHz and 
r+ = 10MHz, as the peak cannot be resolved with the above decoherence rates. 
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FIG. 6: (Colour online) Distribution of system parameters for a 14keV implant strategy where the two donors are implanted 
30nm apart along a 100 crystallographic axis, obtained using a total of 9149 points, a) The distribution of donor separations R. 
b) The distribution of the symmetric energy gap between the ground and first-excited states of the unperturbed qubit, A. c) 
The distribution of the DC field strength required to bring the system into resonance with a AQGHz field, d) The distribution 
of 4>, the angle that parameterises the degree of localisation of the ground-state, e) The distribution of 7^, the resonant Rabi 
frequency. 



